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requirements such as the gauge symmetry, the Lorentz- invariance and 
the mathematical convergence. Therefore, the renormalized results 
(— I ! derived in the subtraction scheme are faithful and have no ambiguity. 

Especially, it is proved that the solution of the renormalization group 
equation satisfied by a renormalized wave function, propagator or ver- 
tex can be fixed by applying the renormalization boundary condition 
' and, thus, an exact S-matrix element can be expressed in the form as 

written in the tree diagram approximation provided that the coupling 
constant and the fermion mass are replaced by their effective ones. In 
the one- loop approximation, the effective coupling constant and the 
effective fermion mass obtained by solving their renormalization group 
equations are given in rigorous and explicit expressions which are suit- 
able in the whole range of distance and exhibit physically reasonable 
asymptotic behaviors. 
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1. Introduction 



In the renormalization of quantum field theories, to extract finite physi- 
cal results from higher order perturbative calculations, a certain subtraction 
scheme is necessary to be used so as to remove the divergences occurring in 
the calculations. There are various subtraction schemes in the literature, such 
as the minimal subtraction (MS)^, the modified minimal subtraction (MS) 
^, the momentum space subtraction (MOM)^, the on-mass-shell subtraction 
(OS) the off-mass shell subtraction (OMS)^ and etc.. However, there 
exists a serious ambiguity problem^'^'^ that different subtraction schemes in 
general give different physical predictions, conflicting the fact that the phys- 
ical observables are independent of the subtraction schemes. Ordinarily, it 
is argued that the ambiguity appears only in finite order perturbative cal- 
culations, while the exact result given by the whole perturbation series is 
scheme- independent. Though only a finite order perturbative calculation is 
able to be done in practice, one still expects to get unambiguous results from 
such a calculation. To solve the ambiguity problem, several prescriptions, 
such as the minimal sensitivity principle^, the effective charge method^ and 
some others^ , were proposed in the past. By the principle of minimal sen- 
sitivity, an additional condition has to be introduced and imposed on the 
result calculated in a finite order perturbative approximation so as to ob- 
tain an optimum approximant which is least sensitive to variations in the 
unphysical parameters. In the effective charge method, the use of a cou- 
pling constant is abandoned. Instead, an effective charge is associated with 
each physical quantity and used to determine the Gell-Mann-Low function in 
the renormalization group equation (RGE). In this way, a renormalization- 
scheme-invariant result can be found from the RGE. In Ref.(9), the authors 
developed a new perturbation approach to renormalizable field theories. This 
approach is based on the observation that if the perturbation series of a quan- 
tity R is not directly computable due to that the expansion coefficients are 
infinite, the unambiguous result of the quantity can be found by solving 
the differential equation = F{R{Q)) where F(R) is well-defined and 

can be expanded as a series of R. Particularly, the finite coefficients of the 
expansion of F(R) are renormalization-scheme-independent and in one-to- 
one correspondence with the ones in the ordinary perturbation series of the 
R. The prescriptions mentioned above are somehow different from the con- 
ventional perturbation theory in which the coupling constant is commonly 
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chosen to be the expansion parameter for a perturbation series. 

In this paper, we wish to deal with the ambiguity problem from a dif- 
ferent angle. It will be shown that the ambiguity problem can be directly 
tackled in the conventional perturbation theory by the renormalization group 
method^"^^^. The advantage of this method is that the anomalous dimension 
in a RGE is well-defined although it is computed from the renormalization 
constant which is divergent in its original definition. In comparison of the 
renormalization group method with the aforementioned approach proposed 
in Ref.(9), we see, both of them are much similar to one another in methodol- 
ogy. This suggests that the renormalization group method is also possible to 
yield the theoretical results which arc free from the ambiguity. The possibil- 
ity relies on how to choose a good subtraction scheme which gives rise to such 
renormalization constants that they lead to unique anomalous dimensions. 
The so-called good subtraction scheme means that it must respect necessary 
physical and mathematical requirements such as the gauge symmetry, the 
Lorentz invariance and the mathematical convergence principle. The neces- 
sity of these requirements is clear. Particularly, in the renormalization group 
approach to the renormalization problem, the renormalization constants ob- 
tained from a subtraction scheme not only serve to subtract the divergences, 
but also are directly used to derive physical results. Apparently, to guarantee 
the calculated results to be able to give faithful theoretical predictions, the 
subtraction procedure necessarily complies with the basic principles estab- 
lished well in physics and mathematics. Otherwise, the subtraction scheme 
should be discarded and thus the scheme-ambiguity will be reduced. Let 
us explain this viewpoint in some detail from the following aspects: (1) 
For a gauge field theory, as one knows, the gauge-invariance is embodied in 
the Ward identity. This identity is a fundamental constraint for the the- 
ory. Therefore, a subtraction scheme, if it is applicable, could not defy this 
identity. As will be demonstrated in latter sections, the Ward identity not 
only establishes exact relations between renormalization constants, but also 
determines the functional structure of renormalization constants; (2) The 
Lorentz-invariance is partly reflected in the energy-momentum conservation 
which holds at every vertex. Since the renormalization point in the renor- 
malization constants will be eventually transformed to the momentum in the 
solutions of RGEs, obviously, in order to get correct functional relations of 
the solutions with the momentum, the energy-momentum conservation could 
not be violated by the subtraction of vertices; (3) Why the convergence 
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principle is required in the renormalization calculation? As we know, the 
renormalization constants are divergent in their original appearance. Such 
divergent quantities are not well-defined mathematically and hence are not 
directly calculable^because we are not allowed to apply any computational 
rule to do a meaningful or unambiguous calculation for this kind of quanti- 
ties. The meaningful calculation can only be done for the regularized form of 
renormalization constants which are derived from corresponding regularized 
Feynman integrals in a subtraction scheme. The necessity of introducing the 
regularization procedure in the quantum field theory may clearly be seen 
from the mathematical viewpoint as illustrated in Appendix A. Therefore, 
in renormalization group calculations, the correct procedure of computing 
anomalous dimensions is starting from the regularized form of renormaliza- 
tion constants. The limit operation taken for the regularization parameter 
should be performed after completing the differentiation with respect to the 
renormalization point. Since the anomalous dimensions are convergent, the 
limit is meaningful and would give definite results. Obviously, the above pro- 
cedure of computing the anomalous dimensions agrees with the convergence 
principle; (4) In comparison with the other subtraction schemes, the MOM 
scheme appears to be more suitable for renormalization group calculations. 
This is because this scheme naturally provides not only a renormalization 
point which is needed for the renormalization group calculation, but also 
a renormalization boundary condition for a renormalized quantity (a wave 
function, a vertex or a propagator), which will be used to fix the solution of 
the RGE for the renormalized quantity. 

Based on the essential points of view stated above, it may be found that 
a renormalized quantity can be unambiguously determined by its RGE and 
thus a renormalized S-matrix element can be given in an unique form with- 
out any ambiguity. To illustrate this point, we limit ourselves in this paper 
to take the QED renormalization as an example to show how the ambiguity 
can be eliminated. As one knows, the QED renormalization has been ex- 
tensively investigated by employing the OS, MS and MOM schemes in the 
previous works^"^'-'^^"^^. But, most of these studies are concentrated on the 
large momentum (short distance) behaviors of some quantities for the sake 
of simplicity of the calculation. In this paper, we restudy the QED renormal- 
ization with the following features: (1) The renormalization is performed in 
a mass-dependent scheme other than in the mass-independent scheme which 
was adopted in many previous works^°~^^. The mass-dependent scheme obvi- 
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ously is more suitable for the case that the mass of a charged fermion can not 
be set to be zero; (2) The subtraction is carried out in such a MOM scheme 
that the renormahzation point is mainly taken to be an arbitrary time-like 
momentum other than a space- like momentum as chosen in the conventional 
MOM scheme. The time-like MOM scheme actually is a generalized mass- 
shell scheme (GMS). The prominent advantage of the GMS scheme is that 
in this scheme the scale of renormahzation point can naturally be connected 
with the scale of momenta and the results obtained can directly be converted 
to the corresponding ones given in the OS-scheme ; (3) The subtraction is 
implemented by fully respecting the necessary physical and mathematical 
principles mentioned before . Therefore, the results obtained are faithful 
and free of ambiguity; (4) The effective coupling constant and the effective 
fermion mass obtained in the one-loop approximation are given exact and 
explicit expressions which were never found in the literature. These expres- 
sions exhibit physically reasonable infrared and ultraviolet behaviors. We 
will pay main attention to the infrared (large distance) behavior because this 
behavior is more sensitive to identify whether a subtraction is suitable or not 
for the QED renormahzation. 

The rest of this paper is arranged as follows. In Sect. 2, we sketch 
the RGE and its solution and show how a S-matrix element can be free 
of ambiguity. In Sect. 3, we briefly discuss the Ward identity and give a 
derivation of the subtraction version of the fermion self-energy in the GMS 
scheme. In Sect. 4, we derive an exact expression of the one-loop effective 
coupling constant and discuss its asymptotic property. In Sect. 5, the same 
thing will be done for the effective fermion mass. The last section serves to 
make some comments and discussions. In Appendix A, we show a couple of 
mathematical examples to help understanding the regularization procedure 
used in the renormahzation group calculations. 

2. Solution to RGE and S-matrix element 

Among different formulations of the RGE (sec the review given in Ref.(14))^° 
we like to employ the approach presented in Rcf.(15). But, we work in a 
mass- dependent renormahzation scheme, therefore, the anomalous dimen- 
sion in the RGE depends not only on the coupling constant, but also on the 
fermion mass. Suppose Fr is a renormalized quantity. In the multiplicative 
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renormalization, it is related to the unrenormalized one F in such a way 

F = ZfFr (2. 1) 

where Zp is the renormahzation constant of F. In GMS scheme, the Zp and 
Fr are all functions of the renormalization point /i = yUge* where is a fixed 
renormalization point corresponding the zero value of the group parameter 
t. Differentiating Eq.(2.1) with respect to the n and noticing that the F is 
independent of /x, we immediately obtain a RGE satisfied by the function F^ 

dF 

/^-^ + 7Fi^ii = (2.2) 

where 7f is the anomalous dimension defined by 

If = n^lnZp (2. 3) 

We first note here that because the renormalization constant is dimensionless, 
the anomalous dimension can only depends on the ratio a — ^p— ^pigji, cr), 
where ttir and g^j are the renormalized fermion mass and coupling constant 
respectively. Next, we note, Eq.(2.2) is suitable for a physical parame- 
ter (mass or coupling constant), a propagator, a vertex, a wave function 
or some other Green function. If the function Fr stands for a renormal- 
ized Green function, vertex or wave function, in general, it not only de- 
pends explicitly on the scale /i, but also on the renormalized coupling con- 
stant Qr, mass rriR and gauge parameter ^r which are all functions of /x, 
Fr = Fr{p, gR{ij,),mR{iJ,),^R{iJ,); ii) where p symbolizes all the momenta. 
Considering that the function Fr is homogeneous in the momentum and 
mass, it may be written, under the scaling transformation of momentum 
p — Xpo , as follows 

Fr{p; gR, ruR, ^) = ^^''^FRipo; gR, ^) (2. 4) 

where Dp is the canonical dimension of F. Since the renormalization point 
is a momentum taken to subtract the divergence, we may set fj, = /iqA where 
A = e* which is taken to be the same as in p = poX. Noticing the above 
transformation, the solution of the RGE in Eq.(2.2) can be expressed as^^ 

FR{p;gR,mR,Ul^o) = A^-e/' ^^-^^^(po; 

gRW,mR{X)X-\^R{Xy,lio) (2. 5) 



6 



where 5'_r(A), mji{X) and ^r(A) are the running coupling constant, the running 
mass and the running gauge parameter, respectively. The solution written 
above shows the behavior of the function Fr under the scaling of momenta. 

How to determine the function Fr(po; • • • , A*o) on the RHS of Eq.(2.5) 
when the Fr(po,---) stands for a wave function, a propagator or a vertex? 
This question can be unambiguously answered in MOM scheme, but was not 
answered clearly in the literature 2'^'^^. Noticing that the momentum po and 
the renormalization point /iq are fixed, but may be chosen arbitrarily, we 
may, certainly, set pi = fi^. With this choice, by making use of the following 
boundary condition satisfied by a propagator, a vertex or a wave function 

FRipo] Qr, niR, ^R, //) \p2=^2= F^\po] Qr, niR, ^r) (2. 6) 

where the function fI^\p; qr, rriR, or) is of the form of free propagator, bare 
vertex or free wave function and independent of the renormalization point 
(see the examples given in the next section) and considering the homogeneity 
of the function Fr as mentioned in Eq.(2.4), we may write 

A^^F^(po; gR{X),mR{X)X-\ UW^f^o) |pg=;.g= F^^\p; ^«(A), m«(A), ei?(A)) 

(2. 7) 

where the renormalized coupling constant, mass and vertex in the function 
Fr\p, ■•■) become the running ones. With the expression given in Eq.(2.7), 
Eq.(2.5) will finally be written in the form 

Fr{p; gR, rriR, ^r) = e/' ^^^^'^F^ (p; gR{X),mR{X),U>^)) (2. 8) 

For a gauge field theory, it is easy to check that the anomalous dimension in 
Eq.( 2.8) will be cancelled out in S-matrix elements. To show this point more 
specifically, let us take the two-electron scattering taking place in t-channel 
as an example. Considering that a S-matrix element expressed in terms of 
unrenormalized quantities is equal to that represented by the corresponding 
renormalized quantities, the scattering amplitude may be written as 

Sfi = <{p[)TUp[:Pi)u%{pi)tDR,,{k)u^^^^^^^ (2. 9) 

where k — p[ — pi = P2 — P2', u'^ip), ^'r{p' iP) and iDR^i,{k) represent the 
fermion wave function, the proper vertex and the photon propagator respec- 
tively which are all renormalized. The renormalization constants of the wave 
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function, the propagator and the vertex will be designated by \/^, ^3 and 
Zy respectively. The constant Zy is defined as 

Zr = Z:^^Z'l^ (2. 10) 

because the vertex in Eq.(2.9) contains a coupling constant in it. 

According to the formula given in Eq. (2.8) , the renormalized fermion wave 
function, photon propagator and vertex can be represented in the forms as 
shown below. For the fermion wave function, we have 

««(P) = e/'^^-(')4°i(P,^«(A)) (2. 11) 

where 

is the free wave function in which mij(A) is the running mass and 

= (2. 13) 

2 a// 

is the anomalous dimension of fermion wave function. For the renormalized 
photon propagator, we can write 



^^«..(^) = e/i'^^^WzDg) (^) 



where 

- -p^t^M. - (1 - ««(A))^] (2. 14) 

is the free propagator with q;r(A) being the running gauge parameter in it 
and 

73(A) In Z3 (2. 15) 

is the anomalous dimension of the propagator. For the renormalized vertex, 
it reads 

r^(p',p) = e/'x7r(A)rg)M(y^^) (2. 16) 

where 
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T^^^'ip',p)=teni\)r (2. 17) 

is the bare vertex containing the running couphng constant (the electric 
charge) eR{X) in it and 

7r(A) = /x-^lnZr = -//-^lnZ2 - ^/x-^lnZ3 (2. 18) 

dfj, an 2 ail 

is the anomalous dimension of the vertex here the relation in Eq.(2.10) has 
been used. Upon substituting Eqs.(2.11), (2.14) and (2.17) into Eq.(2.9) and 
noticing Eqs.( 2.13), (2.16) and (2.19), we find that the anomalous dimensions 
in the S-matrix element are all cancelled out with each other. As a result, 
we arrive at 

(2. 19) 

This expression clearly shows that the exact S-matrix element of the two- 
electron scattering can be represented in the form as given in the lowest order 
(tree diagram) approximation except that all the physical parameters in the 
matrix elements are replaced by their effective (running) ones. For other 
S-matrix elements, the conclusion is the same. This is because any S-matrix 
element is unexceptionably expressed in terms of a number of wave functions, 
propagators and proper vertices each of which can be represented in the form 
as shown in Eq.(2.8) and the anomalous dimensions in the matrix element 
, as can be easily proved, are all cancelled out eventually. This result and 
the fact that any S-matrix element is independent of the gauge parameter 
( This is the so-called gauge-invariancc of S-matrix which is implied by the 
unitarity of S-matrix elements) indicate that the task of renormalization for 
a gauge field theory is reduced to find the running coupling constant and the 
running mass by their RGEs. These running quantities completely describe 
the effect of higher order perturbative corrections. 



3. Ward Identity 

In QED renormalization, the following Ward identity plays crucial role ^ 

AM(p';p)b'=P = -^ (3. 1) 
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where Af^{p',p) represents the vertex correction which is defined by taking 
out a couphng constant e and S(p) denotes the fermion self-energy. Firstly, 
we show how the above identity determines the subtraction of the fermion 
self-energy J2{p)- According to the Ward identity, the divergence in A^{p',p) 
should be, in the GMS scheme, subtracted at a time-like (Minkowski) renor- 
malization point for the momenta of the external fermion lines, p^ = p'^ = /x^ 
which implies ^ = ^' = fJ^- When /j, = m (the fermion mass), we will come to 
the subtraction in the OS scheme. In the case of // 7^ m, the subtraction is 
defined on a generalized mass shell. At the renormalization point /i, we have 

KiP^P) l/=y=/.= (3- 2) 

Thus, the vertex correction may be represented as 

A,ip',p)^L^, + Alip',p) (3. 3) 

where L is a divergent constant depending on /j, and A^(p',p) is the finite 
correction satisfying the boundary condition 

A^(p',p) \^=^=,^0 (3. 4) 

On inserting Eq.(3.3) into Eq.(3.1) and integrating the both sides of 
Eq.(3.1) over the momentum p^, we get 

E(p) - E(/.) = -{^-i^)L- r dp^A'^^ip^p) (3. 5) 

where the momentum po is chosen to make ^0 = Since the last term on 
the RHS of Eq.(3.5) vanishes when p^ — > p^, we may write 

dp^Al{p,p)^{^-fx)Cip') (3. 6) 

-'Po 

where C(p^) is a convergent function satisfying the following boundary con- 
dition 

C(p2) 1^2=^2= (3. 7) 

which is imphed by Eq.(3.4). Substituting Eq.(3.6) into Eq.(3.5) and setting 

E(/x) = A (3. 8) 
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and 

L^-B (3. 9) 

the self-energy is finally written as 

S(p) = ^+(/.-/.)[S-C(p2)] (3. 10) 

where the constants A and B have absorbed all the divergences appearing in 
the The above derivation shows that the subtraction given in Eq.(3.10) 

is uniquely correct in the GMS scheme as it is compatible with the Ward 
identity. According to the subtraction in Eq.(3.3), the full vertex can be 
written as 

where Zi is the vertex renormalization constant defined as 

Z-^ = 1 + L (3. 12) 

and r^{p',p) is the renormalized vertex represented by 

r>',p)=7M + A?(p',P) (3- 13) 

which satisfies the boundary condition 

Tl^^'ip^p) \^=j^=,= l, (3. 14) 

Based on the subtraction given in Eq.(3.10), the full fermion propagator 
may be renormalized in such a way 

where Z2 is the propagator renormalization constant defined by 

Z^^ ^1-B (3. 16) 

and Sp{p) denotes the renormalized propagator represented as 

Spip) = -J (3- 17) 
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which has a boundary condition as follows 

S§{p) |p^=,^= (3. 18) 

In Eq.(3.17), uir and T1r(p) designate the renormalized mass and the finite 
correction of the self-energy respectively. The renormalized mass is defined 

by 

rriR = Z-^m (3. 19) 

where Z^ is the mass renormalization constant expressed by 

Z-^^l + Z2[Am-' + (1 - i^m-^)B] (3. 20) 

Particularly, from Eqs.(3.9), (3.12) and (3.16). it is clear to see 

Zi - Z2 (3. 21) 

This just is the Ward identity obeyed by the renormalization constants. 

Let us verify whether the Ward identity is fulfilled in the one-loop 
approximation. The Feynman integrals of one-loop diagrams in QED have 
been calculated in the literature by various regularization procedures^"^'^^'^^,. 
In the GMS scheme, the fermion self-energy depicted in Fig. (la), according 
to the dimensional regularization procedure, is regularized in the form 

x(l-x)/)- (4-2£)m+ (1 -0(^-2xi))] -2(1-0 
x(l-.)^} (3.22) 



where £ = 2 — |, 



A(p) = p'^x{x - 1) + m^x (3. 23) 



and M is an arbitrary mass introduced to make the coupling constant e to 
be dimensionless in the space of dimension n. According to the definition 
shown in Eq.(3.8) and noticing 

p'^> = {^>-l^)[p'+l^{^>+l^)]+l^' (3. 24) 
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one can get from Eq.(3.22) 



X [2n[l + - 2)x - s{l - x)] - (3 + ^ - 2s)m] 
-2(l-0(l-x)x^^^ (3.25) 



where 



A(//) = - 1) + m^] (3. 26) 

On substituting Eqs.(3.22) and (3.25) in Eq.(3.10), it is found that 

B = [E(p)-A](p-/x)-My=M 



2//' 



2 



x(l - x) - 2(1 - Ox] + ^^y^m - s)x{x - 1) 



\2 



A(/.) 



+5(1 - - 1) + + e - 2£)a;(a; - 1)] - 

x{l-0{l + e){x-lfx^} (3.27) 

For the diagram of one-loop vertex correction shown in Fig. (lb), ac- 
cording to the definition written in Eq.(3.2), it is not difficuh to obtain, in 
the n-dimensional space, the regularized form of the constant L 

-(1 -Oix'-x- 1)//^ + e{l -Oix- 1)/^' - 4m//[(2 - s) 
x(x - 1) + ^(1 - 0(2 - 3a;) + ^e{l - Oi^ ' 1)] + ^'[2(£ 

-1) - (1 - 6)(1 - Oix - 1)]] + (1 - - 

x{m + nf} (3. 28) 
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With the expressions given in Eqs.(3.27) and (3.28), in the approximation 
of order e^, the renormahzation constants defined in Eq.(3.12) and (3.16) 
will be represented as Zi = 1 — L and Z2 — 1 + B. In the hmit £ — > 0, 
these constants are divergent, being not well-defined. So, to verify the Ward 
identity in Eq.(3.21), it is suitable to see whether their anomalous dimensions 
satisfy the corresponding identity 

71 = 72 (3. 29) 

where 7i(/x) = limfr^Q fi^ln Zi{ij,, s) (i=l,2). For the renormahzation group 
calculations, in practice, the above identity is only necessary to be required. 
Through direct calculation by using the constants in Eqs.(3.27) and (3.28), 
it is easy to prove 

7i = 72 

-{ee - 6(3 + Oct + 12^(7^ + 6(3 + e - 2ea) 



(47r) 



2 



2 1 
xaHn^- + A[2^-{3 + 0(r]^- (3. 30) 

0""^ — 1 (7^ — 1 

where a = ^. This identity guarantees the correctness of the one-loop renor- 
malizations in the GMS scheme. In the zero-mass limit (cr — > 0), the identity 
in Eq.(3.30) reduces to the result given in the MS scheme 

^'^^^^8^ (3. 31) 

This result can directly be derived from such expressions of the constants B 
and L that they are obtained from Eqs.(3.27) and (3.28) by setting m=0. 
It is easy to see that in these expressions, only the terms proportional 
give nonvanishing contributions to the anomalous dimensions. However, in 
the case of m 7^ 0, the terms proportional to e^^ in Eqs.(3.27) and (3.28) 
give different results. In this case, to ensure the identity in Eq.(3.29) to be 
satisfied, the other terms without containing e"^ in Eqs.(3.27) and (3.28) 
must be taken into account. 
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4. Effective Coupling Constant 



The RGE for the renormahzed couphng constant may be immediately 
written out from Eq.(2.2) by setting F — e, 

l^'^^M + le{lJ)(^R{lJ) = (4. 1) 

In the above, the anomalous dimension jeil^) as defined in Eq.(2.3) is now 
determined by the following renormalization constant^ 



Ze = ^ = ^3 ' (4. 2) 

Z2Z2 

where the identity in Eq.(3.21) has been considered. The photon propagator 
renormalization constant is, in the GMS scheme, defined by 

= 1 + n(//2) (4. 3) 

where n(//^) is the scalar function appearing in the photon self-energy tensor 
n^ni^ = {kjzi,—k'^gni,)ll{li'^). In view of Eq.(4.2), we can write 

7e = lim jx^ InZe^ ~ lim /x-^ In Z3 (4. 4) 

£^0 d/i 2 £^0 d/i 

For the one-loop diagram represented in Fig.(lc), according to Eq.(4.3), it 
is easy to derive the regularized form of the constant Z^ by the dimensional 
regularization procedure 

Substituting Eq.(4.5) into Eq.(4.4), it is found that 

7^ = —-{1 + 6(7^+ , In } 4.6 

' 12(7r)2^ VI - 4(72 1- VI -4(72-^ ^ ' 

where a = ^. In this expression, the charge e and the mass m are un- 
renormalized. In the approximation of order e^, they can be replaced by the 
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renormalizcd ones e^, and rriB because in this approximation, as pointed out 
in the previous hterature^"^, the lowest order approximation of the relation 
between the e{m) and the eR{mji) is only necessary to be taken into account. 
Furthermore, when we introduce the scaling variable A for the renormaliza- 
tion point and set ^Uq = rriR (which can always be done since the /xq is fixed, 
but may be chosen at will), we have a = ^ = j. Thus, with the expressions 
of Eq.(4.6), Eq.(4.1) may be rewritten in the form 

= P{X) (4. 7) 

where 

/3(A) = -7e(A)eR(A) 

e|(A) 



127r2 



Fe(A) (4. 8) 



in which 



FM = 1 + f + fjW (4. 9) 



(4. 10) 



Upon substituting Eqs.(4.8)-(4.10) into Eq.(4.7) and then integrating Eq.(4.7) 
by applying the familiar integration formulas, the effective (running) coupling 
constant will be found to be 

MA) = ^^feW 

where q;r(A) = an = q;k(1) and 



G{X) = / ^Fe(A) 



= 2 + v^TT - ^ + (I + ^)^^(A) (4. 12) 
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in which 



_ r-VI^cos-^i ifA<2 

" lv0^^cosh-^|, if A>2 '"^^ 

As mentioned in Sect. 2, the variable A is also the scaling parameter of mo- 
menta, p = Apo and it is convenient to put po^ — A*o^ so as to apply the 
boundary condition. Thus, owing to the choice ii^ = m^, we have = '^r'^ 

2 I 

and A = (^)2. In this case, it is apparent that when A = 1, Eq.(4.11) will 

be reduced to the result given on the mass shell, aR{l) — aR — which is 
identified with that as measured in experiment. 

The behavior of the q;k(A) are exhibited in Figs. (2) and (3). For small 
A, Eq.(4.11) may be approximated by 

aR{X) ^ ^A^ (4. 14) 

It is clear that when A — > 0, the aR{X) tends to zero. This desirable behavior, 
which indicates that at large distance (small momentum) the interacting 
particles decouple, is completely consistent with our knowledge about the 
electromagnetic interaction. For large momentum (small distance), Eq.(4.11) 
will be approximated by 

Stt 

This result was given previously in the mass-independent MS scheme. In 
the latter scheme, the (3— function in Eq.(4.8) is only a function of the 
ei?(A) since Fe{X) — 1 due to m = in this case. But, in general, the 
mass of a charged particle is not zero. Therefore, the result given in the MS 
scheme can only be viewed as an approximation in the large momentum limit 
from the viewpoint of conventional perturbation theory. Fig. (3) shows that 
the q;r(A) increases with the growth of A and tends to infinity when the A 

approaches an extremely large value Aq c^^a — e^" (the Landau pole). If 
the A goes from Aq to infinity, we find, the q;r(A) will become negative and 
tends to zero. This result is unreasonable, conflicting with the physics. The 
unreasonableness indicates that in the region [Aq, oo), the QED perturbation 
theory and even the QED itself is invalid^. 
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5. Effective Fermion Mass 



The RGE for a renormalized fermion mass can be directly read from Eq.(2.2) 

I— — A— 



when we set F — niR. Noticing ji-^ — X-fr, this equation may be written in 



the form 

= -7.(A)m^(A) (5. 1) 

where the anomalous dimension 7m(A), according to the definition in Eq.(2.3), 
can be derived from the renormalization constant represented in Eq.(3.20). 
At one-loop level, by making use of the constants A, B and Z2 which were 
written in Eqs.(3.25), (3.27) and (3.16) respectively, in the approximation of 
order e^, it is not difficult to derive 

7™(A) = lim/.|- In = (5. 2) 



where 



F„(A) = 2eA + 6[3 + 2e-^^^ + |] 

12(1 + OA ,,,,,, 3(1 + , 20 1 



here the relation (J — \ has been used. Inserting Eq.(5.2) into Eq.(5.1) and 
integrating the latter equation, one may obtain 

mR{\) = mRe-*W (5. 4) 

This just is the effective (running) fermion mass where niR = mR{l) which 
is given on the mass-shell and 

SiX) ^ ^ £ fan{X)F^{X) (5. 5) 

In the above, the bare charge appearing in Eq.(5.2) has been replaced by 
the renormalized one and further by the running one shown in Eq.(4.11). If 
the coupling constant in Eq.(5.5) is taken to be the constant defined on the 
mass-shell, the integral over A can be exphcitly calculated. The result is 

S{X)^^[M>^)+^M^)] (5. 6) 
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where 



(^i(A) = 3(l-A){^ + [^--i(l + A)]ln 



1-A^ 



and 



(^2(A) = (l-A){(2A-3)^ln 



1-A^ 



3(1 -A) 
A2 



}-2 



in which 



In 



1-A^ 



2[ln(l + A) -tanh"^A], if A < 1 
2[ln(l + A) -cot/i-U], ifA>l 



(5. 7) 



(5. 8) 



(5. 9) 



As wc sec, the function S{\) and hence the effective mass mii{\) are gauge- 
dependent. The gauge- dependence is displayed in Fig. (4). The figure shows 
that for ^ < 10, the effective masses given in different gauges are almost the 
same and behave as a constant in the region of A < 1, while, in the region 
of A > 1, they all tend to zero with the growth of A. But, the mij(A) given 
in the gauge of ^ 7^ goes to zero more rapidly than the one given in the 
Landau gauge {C, = 0). To be specific, in the following we show the result 
given in the Landau gauge which was regarded as the preferred gauge in the 
literature^'^^. 



In the limit A — > 0, 



mii{\) mRe~~i^ = 1.001744mR 



(5. 10) 



(5. 11) 



This clearly indicates that when A varies from 1 to zero, the mK(A) almost 
keeps unchanged. This result physically is reasonable. Whereas, in the region 
of A > 1, the mji{X) decreases with increase of A and goes to zero near the 
critical point Aq. This behavior suggests that at very high energy, the fermion 
mass may be ignored in the evaluation of S-matrix elements. 



6. Comments and Discussions 

In this paper, the QED renormalization has been restudicd in the GMS 
scheme. The exact and explicit expressions of the one- loop effective coupling 
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constant and fermion mass are obtained in the mass-dependent renormal- 
ization scheme and show reasonable asymptotic behaviors. A key point to 
achieve these results is that the subtraction is performed in the way of re- 
specting the Ward identity, i.e. the gauge symmetry. For comparison, it 
is mentioned that in some previous literature^^'^^, the fermion self-energy is 
represented in such a form 

E(p)=A(p2) + (6- 1) 

If we subtract the divergence in the at the renormalization point = /x^, 
the fermion propagator is still expressed in the form as written in Eqs.(3.15) 
and (3.17); but, the renormalization constants Z2 and are now defined 

by 

Z^' = 1 - Z-' = Z2[l + B{^')] (6. 2) 

The one-loop expressions of the A{iJ,^) and B{ii^) can directly be read from 
Eq.(3.22). Here, we show the one-loop anomalous dimension of fermion prop- 
agator which is given by the above subtraction 

72 = lim/x-^lnZ2 

In comparison of the above 72 with the 71 shown in Eq.(3.30), we see, the 
Ward identity in Eq.(3.29) can not be fulfilled unless in the zero- mass limit 
((7 — > 0) . Since the Ward identity is an essential criterion to identify whether 
a subtraction is correct or not, the subtraction stated above should be ex- 
cluded from the mass-dependent renormalization. 

Another point we would like to address is that in the Ward identity 
shown in Eq.(3.1), the momenta p and p' on the fermion lines in the vertex 
are set to be equal. According to the energy-momentum conservation, the 
momentum k on the photon line should be equal to zero, correspondingly, 
the subtraction shown in Eq.(3.2) was carried out at the so-called asym- 
metric points, p'"^ = p'^ = fi'^ and k"^ = {p' — pY = 0. These subtraction 
points coincide with the energy- momentum conservation, i.e., the Lorentz- 
invariance. Nevertheless, the subtraction performed at the symmetric point 
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p/2 _ p2 _ f^2 _ ^g^g Qf^g]^ used in the previous works ^'19,20^ This sub- 
traction not only makes the calculation too complicated, but also violates 
the energy- momentum conservation which holds in the vertex. That is why 
the symmetric point subtraction is beyond our choice. 

As mentioned in Introduction, the GMS subtraction is a kind of MOM 
scheme in which the rcnormalization points arc chosen to be time-like. In 
contrast, in the conventional MOM scheme^'^'^^, the rcnormalization points 
were chosen to be space-like, i.e. = —/i^ which implies — ifj,. In this 
scheme, the one-loop result for the anomalous dimension 7e can be written 
out from Eq.(4.6) by the transformation a — > —ia. That is 

7e = — No il — 6cr H , In , = \ (6.4) 

which is identical to that given in Refs.(6) and (29). In comparison of Eq.(6.4) 
with Eq.(4.6), we see, the coefficients of cr^ in Eq.(6.4) changes a minus sign. 
Substituting Eq.(6.4) into Eq.(4.1) and solving the latter equation, we obtain 
the running coupling constant which is still represented in Eq.(4.11), but the 
function G{X) in Eq.(4.11) is now given by 

2 ^ ^/PT4,2 1 



+v^ln^(l + v^) (6. 5) 

This expression is obviously different from the corresponding one written 
in Eqs.(4.12) and (4.13). At small distance (A oo), Eq.(6.5) still gives 
the approximate expression presented in Eq.(4.15). However, at large dis- 
tance, as shown in Fig. (2), the aR(A) behaves almost a constant. When 
A — > 0, it approaches to a value equal to 0.99986q;r, unlike the ol^I^X) given 
in Eqs.(4.11)-(4.13) which tends to zero. Let us examine the effective mass. 
As we have seen from Sect. 5, the one-loop effective fermion mass given in 
the GMS scheme is real. However, in the space-like momentum subtraction, 
due to ^ = i/x, the effective mass will contain an imaginary part. This re- 
sult can be seen from the function F^(A) whose one-loop expression given in 
the usual MOM scheme can be obtained from Eq.(5.3) by the transforma- 
tion \ i\ and therefore becomes complex. The both of subtractions may 
presumably be suitable for different processes of different physical natures. 
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But, if the effective mass is required to be real, the subtraction at space-hke 
renormalization point should also be ruled out. 

As pointed out in Sect. 4, the effective coupling constant shown in 
Eq.(4.15) which was obtained in the MS scheme is only an approximation 
given in the large momentum limit from the viewpoint of conventional per- 
turbation theory. Why say so? As is well-known, the MS scheme is a mass- 
independent renormalization scheme in which the fermion mass is set to 
vanish in the process of subtraction. The reasonability of this scheme was 
argued as foUows^'^^ .The fermion propagator can be expanded as a series 

1 111111 , , 

- -|- —rm— + — Tin— I'm— -\- ■ ■ ■ (6. 6) 



p — m 

According to this expansion, the massive propagator may be replaced 
by the massless one ^. At the same time, the fermion mass, as the coupling 
constant, can also be treated as an expansion parameter for the perturba- 
tion series. Nevertheless, in the mass-dependent renormalization as shown 
in this paper, the massive fermion propagator is employed in the calculation 
and only the coupling constant is taken to be the expansion parameter of 
the perturbation series. Thus, in order to get the perturbative result of a 
given order of the coupling constant in the mass- dependent renormalization, 
according to Eq.(6.6), one has to compute an infinite number of terms in the 
MS scheme. If only the first term in Eq.(6.6) is considered in the MS scheme, 
the result derived in the this scheme, comparing to the corresponding one 
obtained in the mass-dependent renormalization, can only be viewed as an 
approximation given in the large momentum limit . Even if in this limit, a 
good renormalization scheme should still be required to eliminate the am- 
biguity and give an unique result. To this end, we may ask whether there 
should exist the difference between the MS scheme and the MS scheme ^? 
As one knows, when the dimensional regularization is employed in the mass- 
independent renormalization, the MS scheme only subtracts the divergent 
term having the e-pole in a Feynman integral which is given in the limit 
£ — > and uses this term to define the renormalization constant. While, the 
MS scheme is designed to include the unphysical terms 7 — In 47r (here 7 is 
the Euler constant) in the definition of the renormalization constant. The un- 
physical terms arise from a special analytical continuation of the space-time 
dimension from n to 4. When the two different reno^'malization constants 
mentioned above are inserted into the relation 6 = ^3 ^e^^, one would derive 
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a relation between the two different renormalized coupling constants given 
in the MS and MS schemes if the higher order terms containing the e-pole 
are ignored. It would be pointed out that the above procedure of leading 
to the difference between the MS and MS schemes is not appropriate be- 
cause the procedure is based on direct usage of the divergent form of the 
renormalization constants. As emphasized in the Introduction, according to 
the convergence principle, it is permissible to use such renormalization con- 
stants to do a meaningful calculation. The correct procedure of deriving a 
renormalized quantity is to solve its RGE whose solution is uniquely deter- 
mined by the anomalous dimension (other than the renormalization constant 
itself) and boundary condition. In computing the anomalous dimension, the 
rigorous procedure is to start from the regularized form of the renormaliza- 
tion constant. In the regularized form, it is unnecessary and even impossible 
to divide a renormalization constant into a divergent part and a convergent 
part. Since the anomalous dimension is a convergent function of e due to that 
the factor - disappears in it, the limit e — > taken after the differentiation 
with respect to the renormalization point would give an unambiguous re- 
sult. Especially, the unphysical factor (47r)^ r(l-|-e) appearing in Eqs.(3.25), 
(3.27) and (4.5) straightforwardly approaches 1 in the limit. Therefore, the 
unphysical terms 7 — In 47r could not enter the anomalous dimension and the 
effective coupling constant. Even if we work in the zero-mass limit or in the 
large momentum regime, we have an only way to obtain the effective cou- 
pling constant as shown in Eq.(4.15) in the one-loop approximation. That is 
to say, it is impossible, in this case, to result in the difference between the 
MS and MS schemes and also the difference between the MOM and the MS 
schemes. 

The above discussions suggest that the ambiguity arising from different 
renormalization prescriptions may be eliminated by the necessary physical 
and mathematical requirements as well as the boundary conditions. It is ex- 
pected that the illustration given in this paper for the QED one-loop renor- 
malization performed in a mass-dependent scheme would provide a clue on 
how to do the QED multi-loop renormalization and how to give an improved 
result for the QCD renormalization. 
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Appendix: Illustration of the regularization procedure by a 
couple of mathematical examples 



We believe that if a quantum field theory is built up on the faithful basis 
of physical principles and really describes the physics, a S-matrix element 
computed from such a theory is definite to be convergent even though there 
occur divergences in the perturbation series of the matrix element. The 
occurrence of divergences in the perturbation series, in general, is not to 
be a serious problem in mathematics. But, to compute such a series, it is 
necessary to employ an appropriate regularization procedure. For example, 
for the following convergent integral 

roo 

f(a) = / dxe-""" (A. 1) 

^0 

which equals to ^, if we evaluate it by utilizing the series expansion of the 
exponential function 

^-ax ^ ^ L^^" (A. 2) 

as we see, the integral of every term in the series is divergent. In this case, 
interchange of the integration with the summation actually is not permissible. 
When every integral in the series is regularized, the interchange is permitted 
and all integrals in the series become calculable. Thus, the correct procedure 
of evaluating the integral by using the series expansion is as follows 



f(a) = hm y ^—^ / dxx" 



= lim V ^ — r = lim - 1 

^-*oo,!jr^ n\ n+1 A^oo a 

_ 1 

a 



e 



-aA\ 



(A 3) 



This example is somewhat analogous to the perturbation series in the quan- 
tum field theory and suggests how to do the calculation of the series with 
the help of a regularization procedure. Unfortunately, in practice, we are not 
able to compute all the terms in the perturbation series. In this situation, we 
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can only expect to get desired physical results from a finite order perturba- 
tive calculation. How to do it? To show the procedure of such a calculation, 
let us look at another mathematical example. The following integral 

F{a) = / dx-—-^ (A. 4) 

Jo [x + ay 

where a > is obviously convergent. When the exponential function is 
expanded as the Taylor series, the integral will be expressed as 

F(a) = dxy - — ^{x + aY = / -. ^- / 

Jo ^0 n\ ^ ' Jo {x + ay Jo x + a 

+ dx^^-^{x + ay-^ (A. 5) 

"^0 n=2 ^■ 

Clearly, in the above expansion, the first term is convergent, similar to the 
tree-approximate term in the perturbation theory of a quantum field the- 
ory, the second term is logarithmically divergent, analogous to the one-loop- 
approximate term and the other terms amount to the higher order corrections 
in the perturbation theory which are all divergent. To calculate the integral 
in Eq.(A.4), it is convenient at first to evaluate its derivative, 

dF(a)_ r,,> + ° + 2),-(,.a) (A. 6) 



da Jo {x + ay 

which is cqTial to — e~"/a^ as is easily seen from integrating it over x by part. 
In order to get this result from the series expansion, we have to employ a 
regularization procedure. Let us define 

Then, corresponding to Eq.(A.5), we can write 

dF^(a) „ dx dx ^ (-ly ^ /-A 



- -2 / T^^+l - . V^^(n-2) / dx(x + ay-^ 

da Jo {x + ay Jo {x + ay n\ Jo 



oo (_-\\n -(A+a) -a 

to n\ ' {A + ay a2 
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From the above result, it follows that 



— ^ = hm — = A. 9 

da A^oo da 

This is the differential equation satisfied by the function F(a). Its solution 
can be expressed as 

F{a) = F{ao) - T da^ (A. 10) 

where oq > is a fixed number which should be determined by the boundary 
condition of the equation (A. 9). Now, let us focus our attention on the second 
integral in the first equality of Eq.(A.8). This integral is convergent in the 
limit A ^ oo and can be written as 

dFi{a) , 1 1 1 

^ - hm / dx- = / dx- = - (A. 11) 



da A-^^ooio {x + ay Jo {x + a)^ a 
Integrating the above equation over a, we get 

Fi(a) = Fi(ao)+ln- (A. 12) 

If setting 

Fi(ao) = ln^ (A. 13) 

A* 

where // is a finite number, Eq.(A.12) becomes 

Fi(a)=ln- (A. 14) 

This result is finite as long as the parameter is not taken to be zero and can 
be regarded as the contribution of the divergent integral Fi(a) appearing in 
the second term of Eq.(A.5) to the convergent integral F(a). The procedure 
described above for evaluating the function F(a) much resembles the renor- 
malization group method and the approach proposed in Ref.(9). It shows us 
how to calculate a finite quantity from its series expansion which contains 
divergent integrals. 
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A.l FIGURE CAPTIONS 

Fig.(l) The one-loop diagrams. 

Fig. (2) The one-loop effective coupling constants given in the region of 
small momenta. The solid curve represents the result obtained at time-like 
subtraction point. The dashed curve represents the one given at space-like 
subtraction point. 
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Fig. (3) The one-loop effective coupling constants for large momenta. The 
both curves represent the same ones as in Fig. (2). 

Fig. (4) The one- loop effective electron masses given in the Landau gauge 
and some other gauges. 
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